Bose condensate of photons in dilute nondegenerated atomic gas 
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We considered peculiarities of Bose-Einstein condensation of photons that are in thermodynamic 
equilibrium with atoms of diluted gases. General equations of the thermodynamic equilibrium of 
the system under study were obtained. Solutions of these equations in the case of high tempera- 
tures, when the atomic components of the system can be considered as nondegenerated ideal gases 
of atoms, and photonic component can form a state with Bose condensate, were studied. Transcen- 
dental equation for transition temperature and expression for the density of condensed photons in 
considered system were derived. Analytical solutions of the equation for the critical temperature 
in a number of particular cases were considered. The existence of two regimes of Bose condensa- 
tion of photons, which differ significantly in nature of transition temperature dependence on the 
total density of pumped into the system photons, was revealed. In one case, this dependence is a 
traditional fractional-power law, and in another one it is the logarithmic law. Applying numerical 
methods, boundaries of existence and condensate implementation conditions were determined for 
both regimes, depending on the physical parameters of the system under study. We also showed, 
that for a large range of physical systems that are in equilibrium with photons (from ultracold gases 
of alkali metals to certain types of ideal plasma), the condensation of photons should occur according 
to the logarithmic regime. 
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I. INTRODUCTION 

The phenomenon of Bosc-Einstcin condensation 
(BEC) became nowadays so accustomed that the pos- 
sibilities of its wide applications are seriously discussed. 
It had become a reality due to implementation of BEC 
in various many-particle systems, especially in ultracold 
gases of atoms and molecules (see, eg, [1-5]). As is well 
known, the basis of BEC phenomenon and related phe- 
nomenon of superfluidity is the property of a macroscopic 
number of bosons to occupy the same quantum state. It 
is also well known that such a property of bosons leads 
to a specific form of the distribution function. It was 
first noted by A. Einstein [6] after reading the work of 
S. Bose [7]. Recall that in mentioned paper S. Bose was 
the first to derive an expression for the photon distri- 
bution function. That is, the energy distribution func- 
tion of bosons was originally derived for photons, but 
not for Bose atoms. For this reason, the idea of obser- 
vation of photon BEC seems quite obvious. However, 
until recently one considered to be impossible to carry 
out experimentally the conditions for such condensation 
to be achieved. The reasons for such a statement seem to 
have been caused by several fairly obvious circumstances. 
In the case of a system consisting of Bose atoms, the 
possibility of BEC phase appearance with temperature 
decreasing can be shown even in the ideal gas model. 
In this case, Bose condensate is formed by atoms with 



zero momentum. However it is known that in vacuum 
photons with zero momentum do not exist. It reflects 
currently accepted fact that photon is considered to be 
massless in vacuum. Besides, it is difficult to imagine the 
method to lower the temperature in a weakly nonidcal 
gas of photons. If one considers photons in a medium 
with decreasing temperature, the number of photons is 
reduced due to absorption by substance. But one of the 
BEC transition conditions is the conservation of the total 
number of particles, that is, the non-condensate particles 
and those ones forming Bose condensate. 

It should be noted that there were a few works of the 
last century (see, eg, [8-10]), which in some way touched 
upon the issues related to Bose condensation of photons. 
These works had risen despite of so obvious (as it would 
seem) arguments in favour of photons condensation im- 
possibility. In [8] authors hypothesized the existence of 
a photon Bose condensate in early stages of the Universe 
expansion, assuming a non-zero photon mass. Authors of 
[8] also considered the possibility that condensate evapo- 
ration due to interaction with charged particles leads to 
the formation of longitudinal relict radiation. It was as- 
sumed that this radiation forms the main part of Universe 
mass. In [9] authors studied the process of establishing 
equilibrium between radiation and matter in fully ion- 
ized plasma. Solution of kinetic equation indicated the 
possibility of photon BEC formation in the system with- 
out photon absorption. Both the scattering processes 
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(for example, photons by electrons) and the processes of 
emission and absorption were considered to be respon- 
sible for the relaxation mechanism. However, it was ar- 
gued in [9] that the presence of photon absorption leads 
to the replacement of photons condensation phenomenon 
to its temporary accumulation at the region of low fre- 
quencies. In [10] authors studied both the possibility of 
formation of effectively two-dimensional photonic super- 
fluid in Fabry-Perot cavity, and the existence of Bogoli- 
ubov elementary excitations in the system under study, 
as well as the methods of experimental observation of 
such superfluidity, based on Landau criterion. 

Until recently such works were considered to be rather 
"exotic". The situation had dramatically changed after 
appearance of the article [11] (see also [12]), in which the 
implementation of BEC of photons in a real experiment 
at room temperature was announced. In this work, the 
photons are not allowed to leave the system (that is a spe- 
cial dye) with the help of specifically placed mirrors with 
extremely high reflectivity. The system was pumped by 
photons with the help of external laser. Relaxation mech- 
anism is determined by multiple absorption of laser pho- 
tons by dye molecules and re-emission of a lower energy 
photons. Thus it became possible to achieve the equi- 
librium state of photons and environment at room tem- 
perature without photons loss. Further pumping photons 
into the system led to the increase of the density of " free" 
photons that created the possibility of its condensation. 
Note that in mentioned paper the system under study is 
also considered to be effectively two-dimensional. As for 
the mass of the photon, the calculations of authors [11] 
gives numerical value approximately equal 6.7 • lCP 33 .g. 

It is effective mass of photons we are talking about. It 
appears due to a significant change of the photon disper- 
sion law in matter as compared with its linear dispersion 
law in vacuum. As is well known, the spectrum of pho- 
tons in matter reveals so-called 'cut-off frequency', that 
is the final value of the frequency at zero wave vector. 
The existence of a finite cut-off frequency is analogous to 
the presence of an effective mass in the photon dispersion 
law (see further). We emphasize that the appearance of 
cut-off frequency in spectrum and, therefore, the effective 
mass of photon, is caused by its interaction with matter. 
In this sense photon in matter can be considered as a 
quasiparticle. Note also that photon effective mass can 
appear as a result of boundary conditions. For example, 
if the system is situated between two concave mirrors 
with high reflectivity, as it is in case of real experiment 
[11,12]. The presence of these mirrors does not only con- 
serve the total number of photons in system and make 
it effectively two-dimensional, but also contributes to the 
establishment of standing wave along a perpendicular to 
the mirrors' surface. The latter circumstance leads to 
the appearance of frequency cut-off in two-dimensional 
photonic spectrum (see in this regard [10]). 

Thus, the conditions of the experiment [11] have elim- 
inated all mentioned above difficulties , that were earlier 
preventing to realize the photon BEC experimentally and 



to implement the "scientific sensation". 

However, an achieved breakthrough enables, in our 
opinion, to go beyond the framework of photon BEC 
studies in effectively two-dimensional systems, and to 
look for the models of other physical systems, where 
the implementation of photon Bosc condensation is also 
available in the region of relatively high temperatures. 
In the present paper, in this regard, we have deter- 
mined conditions for Bose-Einstcin condensation of pho- 
tons which are in thermodynamic equilibrium with the 
atoms of diluted gases. Special attention is paid for the 
case of high temperatures, when rarefied atomic gases 
can be regarded as nondegeneratc. 



II. GENERAL EQUATIONS OF 
THERMODYNAMIC EQUILIBRIUM OF 
PHOTONS AND ATOMS OF DILUTED GAS 

To solve this problem it is proposed to proceed from 
the following model. Consider the system to be composed 
of an ideal gas of Fermi or Bose atoms that are in ther- 
modynamical equilibrium with photons. Atoms of this 
gas can be only in two quantum states: the ground one 
and the excited one (two- level atoms). In other words, 
the transition between one atomic state to another one 
happens as a result of absorption or emission of a pho- 
ton. Thus, the excited atom can be regarded as a bound 
state of a photon and an atom in the ground state (in 
this connection see [13]). The possibility of considering 
the system as an ideal gas of atoms and photons is based 
on the assumption that once thermodynamic equilibrium 
is reached in a system, the interaction between its com- 
ponents can be neglected. 

In accordance with mentioned above, we actually put 
into consideration two types of atoms residing in equilib- 
rium with a gas of photons. We assume that type "1" of 
atoms is characterized by a set of quantum numbers «i 
and the distribution function of atoms 



*i (p) = jexp 



(p) - Mi 
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(e Ql is the energy of levels of an atom at rest with a set 
of quantum numbers a\, e ai < 0, [i\ is the chemical 
potential, m is the mass of atom) and the type "2" of 
atoms is characterized by a set of quantum numbers 
and the distribution function of atoms 
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e a2 (p) = e a2 + 



2m' 



where e a2 is the energy of levels of an atom at rest with 
a set of quantum numbers a2, e Q2 < 0, and fi2 is the 
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chemical potential. In formulae ([lj - pi), the sign "— " 
must be chosen in the case of bosonic gas of atoms, and 
the sign "+" in the case of fermionic gas. We consider 
also photons that are in thermodynamic equilibrium with 
two-component diluted gas, to be described by the dis- 
tribution function 



fph (k) = jexp 



huj (k) 



T 



- 1 



(3) 



where lj (fc) is the dispersion law of photons and /i* is 
their chemical potential. The introduction of non-zero 
chemical potential of photons implies the conservation 
of their total number in a given thermodynamic state of 
system. 

Condition of thermodynamic equilibrium means that 
at given temperature in system, the atomic N and the 
photonic N p h total numbers are conserved: 

N = 9ai ( p ) + f« 2 /a a (P) > (4) 

P V 



/ill 



= 9* Y fp h + Y f<** ( p ) • 



(5) 



However, as is easily seen, before solving the problem 
it is necessary to clarify the specific type of photons dis- 
persion law cj(fc). As we mentioned above, the type of 
dispersion law determines the effective mass of photons 
in matter, thus determining the conditions of their Bose 
condensation. In the present paper, following [9-12], we 
assume photon dispersion law to be quadratic on the 
wave vector (or photon momentum). It can be 'justi- 
fied' if one starts from the expression for energy of the 
photon in matter as a relativistic object [10, 12]: 



hiu (k) 



lj 2 + v 2 k 2 , 



(7) 



where v is the speed of light in matter and ujo is the 
spectrum cut-off frequency. Then in the region of small 
wave vectors one obtains 



hw (k) = h\/ujQ + v 2 k 2 fa huj 1 
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2 



v 2 {hkf 

2Hljq 



= Hujq + 



(hk) 
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which implies that photons with momentum hk possess 
an effective mass m* 



where <? Ql , g a2 are parameters that take into account the 
degeneracy multiplicity of atomic levels with the set of 
quantum numbers ai, ui respectively. For example, if 
the energy level with a set of quantum numbers ol\ is 
spin-degenerated, then g ai = 2S ai + 1. Similarly, the 
quantity g* takes into account the possible energy de- 
generacy of a photon with energy hu (k). For example, 
if this degeneracy is associated only with photons polar- 
ization, the sum over polarizations gives g* — 3. Recall 
that in the presence of matter photons can have a lon- 
gitudinal polarization, and it was taken into account in 
the previous statement. We emphasize that equation ([5| 
reflects the conservation of the total number of pumped 
into system photons, that is, the total amount of free 
photons and photons absorbed by two-level atoms. Since 
atomic states ai, ai differs in single-photon transition, 
then chemical potentials /i l5 /i 2 and fi* are supposed to 
be connected by following equation (so-called Gibbs rule 
or chemical reaction condition, see in this regard [13] and 
[14]): 



Hl+H = [12- 



(6) 



Equations Q - (6]) form a complete system of cou- 
pled equations to describe thermodynamic equilibrium 
between photons and ideal gas of two-level bosons or 
fermions. Deriving from these equations chemical poten- 
tials as functions of temperature, total number of photons 
and total number of atoms, thus we determine in accor- 
dance with ([lj - ([3]) the distribution functions of all the 
three components of the system under study. 
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As was already mentioned, the presence of cut-off fre- 
quency uiq in photon spectrum may be caused by differ- 
ent reasons, both related to the interaction of photons 
with matter, and the presence of external constraints 
such as mirrors. In this sense, the shape of spectrum 
([7]) is not dictated by "relativism" in any way. In fact, as 
is well known, the dispersion equations for electromag- 
netic waves in plasma yield the following dispersion law 
for transverse waves (see e.g. [15]): 



u? t (k) =ojq+ c 2 k 2 , 



(10) 



where c is the speed of light in vacuum, and ojq is the 
plasma frequency (the frequency of Langmuir oscilla- 
tions) 



4nne 2 



(11) 



In the last formula, n is the density of electrons in plasma 
and e is the elementary charge. Note that longitudinal 
waves in plasma also have a quadratic dispersion law 
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uji (k) = cj ( 1 + - (kr D ) 



(12) 



where Td is the radius of Coulomb screening (Debye ra- 
dius) 
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(13) 



One can see that in the case of fully ionized plasma the 
existence of cut-off frequency is not caused by object's 
relativistic nature at all. 

Naturally, the consistent description requires the calcu- 
lation of cut-off frequency for the photons' dispersion law 
and determining the speed of propagation of electromag- 
netic waves in the system under study. However, both the 
derivation of dispersion equations for specific mediums, 
and their solutions are separated problems, and solving 
each of them involves considerable mathematical difficul- 
ties. Such difficulties, that rise while deriving and solving 
the dispersion equations for electromagnetic waves in the 
systems with bound states and in the presence of a Bose 
condensate, can be estimated from works [16-19]. There- 
fore, in this paper we shall not intend to determine the 
cut-off frequency of the photon spectrum and their veloc- 
ities in the present system, considering them to be given. 
Moreover, in all subsequent calculations, the energy of 
the photons (see (l7|)) to be written as: 



huj (k) = fku (p) = Hujq 



P 

2m* 



(14) 



where the notation for the photon momentum p = hk 
was introduced and the quantity m* is given by (9 1. 

Taking into account (14) and proceeding in (4) and 
([5]) from the momentum sums to integrals according to 
the standard rules of quantum mechanics, we obtain the 
following general equations of thermodynamic balance of 
radiation and matter (ideal gas of atoms): 
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where n and n p h are total densities of atoms or photons, 
respectively: 



N 



Nph 
V 



(18) 



(V is the volume of the system). 

It should be noted that in the framework of our model 
the values of n p h and n should be in the same order 
of magnitude. This follows from the fact that we have 
declared a specific system of two-level atoms. This as- 
sumption is closer to reality if one considers 



Nph ~ N, 



n ph 



(19) 



In this case one can hope that the number of excited 
atoms with the energy levels higher than e Q2 , is negligi- 
ble. 



Equations (15) - (18) should be considered as source 



equationes for solving assigned in the present article 
problem: to determine the conditions under which there 
is a Bose condensate in a gas of photons that are in ther- 
modynamic equilibrium with the atoms of dilute gases. 
Analytical solution of these equations is not generally 
possible. However, the most interesting case is the ap- 
pearance of a Bose condensate of photons in the system 
under study in conditions when the atomic components 
are far from degeneration. In this case the system of 
equations (15) - (17) simplifies greatly, and it is possible 



to obtain some particular solutions in analytical form. 



III. SOLUTION OF THERMODYNAMIC 
BALANCE EQUATIONS FOR THE SYSTEM IN 
SYMMETRIC PHASE 



First of all, we study the solutions of equations (15) 



(17) in the symmetric phase of the system far away from 



the transition point, when all the three components of 
the system can be considered as remote from degenera- 
tion. It will help us to determine relations between the 
system parameters that contribute to achieving the max- 
imum possible number of free photons with temperature 
decreasing. Such behaviour of the photon component 
would maximize the transition temperature of photons 
to a state with BEC. Recall that in the theory of phase 
transitions it is accepted to call the phase of higher sym- 
metry just as "symmetric phase", and the phase of lower 
symmetry as " asymmetric phase" . 

To begin, we consider the case where all the three com- 
ponents are far from degeneration. This case is interest- 
ing due to the following circumstance: we are looking for 
the situation when in nondegenerated state the density of 
free photons increases with decreasing the temperature. 
This behaviour of the photon component would effec- 
tively help in experimental realization of photonic BEC 
in the system under consideration. Conditions of non- 
degeneracy of all the three components can be written 
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cxp I I > 1, exp I ) > 1, 

/ ?VJJQ — /i* 
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(20) 



Then momentum integrals from (151, (16) can be calcu- 
lated 
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and equations of thermodynamic equilibrium (151 - (17) 
can be written as: 



n = e^/ T g ai (T)+e^^ T g a2 (T) , 
n ph = e»'/ T g* (T) + e^/ T g a2 (T) , 

e t*i/T e n*/T = e M2/T^ 

where the following notations were introduced 



(22) 



n ph {T) = e^ T g* (T) 



(25) 



where we have introduced the tilde sign to avoid a con- 
fusion with the density of total pumped photons n p h. 

Note also that in the case of nondegenerated atomic 
components, their statistics has a weak influence on the 
thermodynamic state of the system, since in the leading 
approximation both bosonic and fermionic distribution 
functions have the same (Boltzmann) form. 

Solution of the system of equations ([22]) can be given 
as: 



n ai (T) = 
n a2 (T) 



[1 + A(T)Y 
A (T) n 



[1 + A(T)Y 
n ph (T) =k(T)A (T) , 



(26) 



where, to simplify the form of expressions, we introduced 
following notations: 



k(T) 



9* (T)g ai (T) 
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(28) 



K 2 (T) K (T) 

Consider the expression for the density of free photons 



(see (26)), taking into account expression (28): 
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cxp(-£ Q2 /T) = g a2 (T) , 
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exp(-e ai /T) = g ai (T) , 
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exp (-huo/T) = g* (T) . 



Note that quantities n ai (T), n a2 (T), 

n ai (T) = e^' T g ai (T) 
n a2 (T) = e^' T g a2 (T) 



(23) 



(24) 



represent the density of atoms in quantum states a,\ and 
«2j respectively, whereas the quantity n p h (T) is the den- 
sity of free photons in the system: 



n P h (T) = - [(n ph -n)-K (T)] 
(n ph -n) 2 + 2 (n ph + n) k (T) + k 2 (T) 



1/2 



(29) 



It is clear that the temperature dependence n P h (T) is 
determined by a single parameter n (T) (see p7|). To 
determine the optimal values of the system parameters 
that would implement the maximum possible number of 
free photons with decreasing the temperature, we shall 
involve numerical calculations. To do it, we introduce 
the dimensionless quantities 



T 



e a . 2 — e ai — hu>o 



i g P hg ai ( m 



n p h 9a 2 V 27r ^ 2 
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n p h 
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in terms of which the expression ( 28 1 for the density of 



free photons can be reduced to the form: 



considered. Due to this circumstance, the inequality (32) 
has to be written correctly in the following form: 



n P h (r) 1 

— = --<,a 

n ph 2 

a 2 + 2(a-2)6>|T| 3/2 exp 



I |3/2 

\t\ 1 exp 



? 2 |T| 3 exp 
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Figure 1 shows 3D plot of the ratio n p h (r, 9) / n p h as a 
function of dimensionless parameters r and 9. This plot 
shows that the density of free photons increases with de- 
creasing the temperature (dimensionless parameter |r|) 
in the range of < r < \. Moreover, this tendency is 
observed almost for all 9. Thus, analyzing the expres- 



sion (31 1 (or (28)), we can conclude that if the physical 
characteristics of the system satisfy the relation 



< 



T 



hujQ ^ 2 ' 



(32) 



then decreasing the temperature the number of free pho- 
tons in system increases, that promotes the appearance of 
photons BEC. Here we must make the following remark. 



The presence of zero in the left-hand side of ( 32 ) is justi 



fied only in a formal sense. As temperature decreases, the 
system approaches to the critical point, where equations 
(22), obtained under the assumptions that all the com- 
ponents of the system to be non-degenerated, T ^> T c 
(see (20)), become incorrect. For this reason, a tempera- 
ture below the critical value in condition ( 22 ) cannot be 




FIG. 1: Dependence of the density of free photons in sys- 
tem (dimensionless quantity r\ = n p h/n p h, see 31 1) on tem- 
perature r and describing parameter 9, see ( 30 ijplotted for 



a = 0.5 (see also (30 1). The density of free photons increases 
with decreasing the temperature in the range < r < 1/2. 
This result does not depend on 6 and depends negligibly on a. 
Therefore, it leads to condition (331 (see also 32), facilitating 
Bose condensation of photons in dilute nondegenerated gases 
of atoms. 



> 1 



2T 



T»T C 



(33) 



Using solutions (126]) and taking into account (27) 



28b, one can justify that in the case when (33) is satis- 



fied, conditions (20) for all the three components to be 



nondegenerated are satisfied too. It is necessary to note 
here the following circumstance. Expression (8) repre- 
sents the closeness of the photon frequency u> (k) to the 
frequency ljq . If thus one assumes that the frequency ujq 
is close to the resonant frequency 



foujQ 



then the implementation of condition (33 1 is only possible 
in the case of <C 1. Note also that this relation 
holds in the real conditions of current experiments on 
the condensation of photons [11, 12], where ~ ~. 



IV. BOSE CONDENSATE OF PHOTONS IN 
DILUTE NONDEGENERATED GAS OF ATOMS 

We now proceed to the study of equilibrium state of 
system with the presence of Bose condensate of photons. 
The most interesting case reveals when the BEC of pho- 
tons appears in conditions of nondegeneration of atomic 
components in both quantum states. In such a situation, 



conditions ( 20 ) should be written as: 
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therefore the first two of the expressions (21) remain 



valid. Then using (17), equations (15) and (16) can be 
reduced: 
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( mT V 

Expressing the quantity exp (/l*2 /T) from the first equa- 



tion of ( 35 ) and substituting it into the second one, we 
obtain tl 



ne following equation 
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ng a2 exp [-^f 2 -] 



(36) 



g ai exp 



T 



9a 2 exp [-^j 



which need to be supplemented by the equation ( 17) con- 



necting the chemical potentials of all the three compo- 
nents. Recall that the left-hand side of the equation (36) 



is an expression for the density of free photons in the 
system. In this sense, the equation (36) is similar to 



the source equation for the study of BEC in an ideal 
gas of atoms, see for example [14]. In other words, we 
have reduced the problem of Bose-Einstein condensation 
of photon gas, which is in thermodynamic equilibrium 
with nondegenerate atomic gas, to the problem of Bose- 
Einstein condensation in an ideal gas of photons with 
a certain effective density. Thus the right-hand side of 
equation ( 36 ) should be considered as such effective den- 
sity. 

As is well known (see e.g. [14]), for an ideal gas of Bose 
atoms with the binding energy e a in the temperature 
range T < T c , the chemical potential \x (T) must satisfy 
the condition /x (T) — e a . The same condition should 
define the transition temperature of BEC in ideal gas, 
/J,(T C ) = e a . Thus, in the case of ideal gas of photons 
according to (36), we have (see also [13]): 



fi* (T < T c ) = fku . 



(37) 



Taking into account (37) and following the procedure of 



[14] (see also [13]), one obtains from (36) an intricate 
transcendental equation for the critical temperature T c 
as a function of the total density of photons n p h and 
atoms n in the system, 



exp 



x < g ai exp 



M2 

T 



mT 
2irh 2 



-3/2 



T 



9a 2 exp 



T 



_! (40) 



and, taking into account the equation (17) rewritten in 
form 



fix + huj a = /i 2 , 

also the chemical potential [i\ 



(41) 



exp 



Ml 
T 



x < g ai exp 



mT 
2nh 2 

hujQ 



-3/2 



T 



exp 



9a 2 exp 



_r (42) 



Thus, in accordance with (22), the densities of atomic 
components are determined. 



In some particular cases, the expression (38) for the 



critical temperature can be significantly simplified. In 
fact, if 



< 



exp 



°ai 



T c 



< 1, 



(43) 



the critical temperature, as is easily seen, is given by 
formula: 



2/3 



2m* 



<?*r(3/2)C(3/2) j ' 



(44) 



g* (2m*T c ) 3/2 



4tt 2 H 3 
n p h — n ^ 1 H exp 

9 a 2 



r (3/2) C (3/2) 



-■OL2 °ai 



hujQ 



-i (38) 



and the expression for the density of condensed photons 
n P °h d {T) as a function of temperature T < T c , for the 
same values of n p u and n: 



n c p ° h nd (T) = n ph -n\l + ^- exp 

9a 2 



-^r (3/2) C (3/2) 



-a 2 °ai 



hujQ 



T 

3/2 



(39) 



Expression (37) yields both the chemical potential ^ 



from the equation ( 35 ) 



"ph 



9ct± + 9a 2 



(45) 



As is easily seen, the expression for the critical temper- 
ature is identical by its form to the expression for the 
temperature of Bose-Einstein condensation in an ideal 
gas of ordinary Bose atoms with the mass m* and the 
density equal to n e JJ . In this limiting case, the transi- 



tion temperature (44) depends on the density of atomic 



and photonic components, as well as on degeneracy mul- 
tiplicity of the atomic levels. The obtained expression 



(44) for the temperature T c must satisfy the conditions 
(43) 



< 2m* (e a2 — e ai — hioo) 



g T(3/2)C (3/2) ] 
4ir 2 h 3 n e p { f J 



2/3 



< 1. 



(46) 



This inequality links all the physical characteristics of the 
system under study. 

In the case when parameters of structural units of the 
system satisfy relations 



-Ct2 c Qi 
~fr. 



< exp 



< -1, 

a 2 ~~ £ a x — 



This expression has the same appearance as if photon 
Bose condensation occurs in an ideal gas of photons in 
the total absence of atoms in the system. However, it 
is necessary to keep in mind that the relation n p h ~ n 
(see (19)) has to be satisfied always , since we are in the 
two-level atoms framework. Nevertheless, assuming that 
in all considered above limiting cases (see (43) - (51)), 



(47) 



< 1, 



the density n p h is the same, then in the last case (49), 
the critical temperature, given by expression (51), will 
be the highest. 



eff 
ph 



Starting from (40)-(42) one can also verify that the 



the critical temperature is also given by the formula (|44|, relations p4|, which determine the conditions of atomic 
but with other expression for n, 



eff _ 
1 ph ~ 



^ph 



(48) 



It is easy to see that in this limiting case, the dependence 
on atomic levels degeneracy multiplicity vanishes. Fur- 
thermore, considering the necessity for condition n p h ~ n 
to be satisfied (see (19)), it is easy to conclude that the 
transition temperature in this case in accordance with 
(38), (48) might be extremely low. However, it could be 



guessed also by examining Fig.l for negative values of r. 
In the case of correctness of the inequalities 



components to be nondegenerated, are always satisfied. 
The reason, ensuring their implementation, is the exis- 
tence of a strong inequality m* -C m. 

Expressions (38) and (39) allow solving the inverse 



problem, i.e. what should be the relationship between the 
physical parameters of the system for photon condensa- 
tion to occur for given, , and desirably " high" , tempera- 
tures? The term "high temperature" in the present arti- 
cle assumes the temperature when atomic components 
are degenerated, including, for example, the range of 
room temperatures. For this purpose, we consider again 
the expression ( 38 ) , rearranging it to the form: 



huin 



T, 



> 1, 



exp 



Hujq 



Tr 



(49) 



» 1, 



the critical temperature is also given by theformula ( 44 1 , 



and one can easily convince that in this case effective 
density equals 



eff _ 

% ph " 



(50) 



In this limiting case the dependences on the atomic com- 
ponent density, as well as on atomic degeneracy multi- 
plicity, vanish. Again, in the last two cases it is also nec- 
essary for conditions (47) and (49) to be satisfied. In this 



case we get the relationship between the system physical 
parameters, similar to (46). 



It is necessary to resurrect here the question of the 
optimal relation between the parameters of the system 
structural units. We have shown previously (see (32) 



and (33)), that if T ^> T c , the most advantageous is the 
condition (e a2 — e ai — huo) > 2T, when with lowing the 
temperature the number of free photons in the system 
increases, that promotes the appearance of photons BEC 
at higher temperatures. The last case allows us to verify 
it directly. In fact, in this case the transition temperature 
in accordance with (44) and (50) is given by expression: 



4:TT 2 n ph 



2m* {g* T (3/2) C (3/2) 



2/3 



(51) 



+n< 1 



9* (2m*T c 
4tt 2 ^ 3 



\3/2 



goi 
9a 2 



exp 



T (3/2) C (3/2) 



hujQ 



Tr 



(52) 



As an example, we examine now the solutions of an- 
nounced problem in two (may be said extremely differ- 
ent) cases. In the first case, we consider conditions of 
appearance of Bose condensate of photons that are in 
the thermodynamic equilibrium with gases, which have 
parameters that are close to those of ultracold vapors of 
alkali metals. In the second case, we consider the same 
issue for photons that are in equilibrium with "classic" 
ideal gas of atoms at room temperature. 

To begin, we examine the possibility of appearance of 
photon BEC in the first case. To do it, one has to assign 
the values of density and temperature in (52). Choosing 



these parameters, it is natural to focus on their values in 
real experiments of trapped BEC in alkali metal vapors 
[1-3]. For this reason, we shall assume n ~ 10 13 cm -3 
and T ~ 10~ 5 i^. For such temperatures and densities, 
the vapors of alkali metals should remain nondegenerate. 
One can estimate the first term on the right-hand side of 
(52), assuming m* ~ 10~ 33 .g (as in [11]) 



9* (2m*T c 
4tt 2 H 3 



,3/2 



T (3/2) ( (3/2) - 1 cm- 



Estimating the second term, we can use one of the basic 
requirements of the present article, namely, the interrela- 
tion n p h ~ n to be correct. Then from (52) one obtains: 
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n p h ~ 1 cm + 



10 13 cto -3 



1 + ex P t( £ "2 ~ £ «i - &*>o) /^ c ] 



Thus, to satisfy the condition n p h ~ n, it is necessary to 
satisfy the following relation: 



(53) 



Recall that the condition n p h ~ ra was introduced in or- 
der to justify the use of the assumption of two-level atoms 
(see (19)). As it follows from these estimates, the contri- 
bution of the first term in (52 1 can be neglected. Thus, 



one concludes that the Bose condensation of photons in 
ultracold vapors of alkali metals can be possible if the 
following approximate equation is held: 



n p h 



1 + exp 

9a 2 



huio 



(54) 



We now consider the possibility of Bose condensation 
of photon gas that is in thermodynamic equilibrium with 
an ideal gas of atoms at room temperature, T ~ 300K. 
We assume that the density of the atomic subsystem is 
close to the density of the "classical" ideal gas, that is n ~ 
10 19 cm -3 . Again, for these conditions one can estimate 
the first term in ( 52 1 as 



9* (2m*T c 



,3/2 



-r (3/2) C (3/2) - 10 u cm" 



4tt 2 /i 3 

As in the previous case, estimating the second term we 
shall keep in mind the condition n p h ~ n. Then from 
(52 1 one gets: 



n p h 



10 11 cm' 3 



1 + exp [(e c 



huj ) /T c 



Thus, as in the previous case, it turns out that the re- 
quirement n p h ~ n leads to the relation ( 53 ) . Further- 



more, as in the previous case, the contribution of the first 
term in (52 1 can be neglected, and again one arrives to 



the expression ( 54 ) 



The result is really amazing, taking into account the 
fact how significantly the systems under consideration in 
two cases differ. Of course, it should be understood that 
the result of the above calculations depends on the value 
of the photon effective mass m* , which in various systems 
may differ significantly from the value 6.7 • 10 _33 <?, given 
in [11]. However, we should emphasize here that the 
value to* ~ 10~ 33 g, taken for above numerical estima- 
tions, corresponds to the most interesting case of the vis- 
ible light, i.e. wavelength range A € {380 nm; 750 nm}. 
As was mentioned above, the effective mass of photon 



is determined by the spectrum cut-off frequency uiq and 
the speed of propagation of photons v in the medium (see 
Q). These values also determine the form of photon's 
dispersion law in the medium (see in this connection |7]) 
and ([8|). Therefore, mentioned above estimation results 
clearly indicate the fact that there are conditions un- 
der which one can neglect the contribution of the "tradi- 
tional" for ideal gas term in the expression for the critical 
temperature (i.e. the first term on the right-hand side 
of(52)). Consequently, there should be physical systems, 



similar to those ones considered before, where expression 
for the transition temperature differs cardinally from the 
" traditional" dependence for the case of ideal gas of or- 
dinary atoms. In the next section, we shall consider such 
situation in details. 



V. MECHANISM OF FAST CONDENSATION 
OF PHOTONS IN DILUTE NONDEGENERATE 
ATOMIC GAS 



To examine the issue that had risen in the previous sec- 
tion, let us return to the equation (52 1 being rearranged 
as follows: 



^ =T(3/2)C(3/2) 

9ai 



+ 1 



Set 2 



exp 



9* (2to* T c 
4:Tt 2 h 3 n 
e a . — huio 



,3/2 



(55) 



As was noted earlier, the first term in ( |55| corresponds 
to the traditional character of Bose condensate behavior 
that occurs gaining the phase transition temperature, de- 
fined, for example, by the expression (51 ) when the con- 



dition (49 1 is true. This character of Bose condensation 



has been observed experimentally as in ultracold gases 
of alkali atoms (see e.g. [1-5]), as well in systems with 
conserved number of photons [11,12]. The second term of 
( 55 1 gives a significantly different behavior of the critical 



temperature of Bose condensation of photons in dilute 
nondegenerate gas of atoms. We shall reveal this fact in 
details. 



We stress once again that ( 55 ) is a transcendental equa- 



tion with regard to the temperature of Bose condensation 
T c , and it can not be solved in elementary functions. On 
this account, in Section 4 we were only interested in the 
various limiting cases, which (in our opinion) are close to 
the current experimental conditions for the photon con- 
densation in the presence of matter [11]. However, since 
the system under study differs significantly from the sys- 
tem of work [11], one should also consider the limiting 
case, when the second term of (551 is dominant. Some 



estimates that indicate the legitimacy of the existence of 
such a regime were given in the end of the previous sec- 
tion. In this limiting case, the expression (55 1 simplifies 
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n p h 
n 



exp 



n ph < n. 



C &1 



IlLUn 



(56) 



We specify here that in this section we consider posi- 
tive values of the energies difference e a2 — e ai — Itwq, see 



in this regard comments to the formula (48) and Fig.l 



Accordingly, the critical temperature of the photon Bosc 
condensate phase formation in the system under study 
for this regime is given by 



9o 



n p h 



9al n ph 



(57) 



Note that the expression (57) reveals the logarithmic 



dependence of Bose condensation temperature in the sys- 
tem under study on the densities of its components. This 
behavior significantly differs from traditional power de- 
pendence (for example, see ( 51 )). In this case, the feature 
of the logarithmic regime of condensation is the following 
circumstance. With a proper pumping the system with 
photons, the critical temperature of Bose condensation 



in accordance with formula (57) can unlimitedly increase 



within a finite density of photons. From a physical point 
of view, it means that in this regime the Bose condensa- 
tion of photons can take place much faster than in men- 
tioned earlier case of the power dependence (51). In this 



regard, it makes sense to talk about the mechanism of 
fast condensation of photons within the conditions of the 
regime (57). As far as we know, previously in literature 



this kind of dependence was not mentioned. 

One can easily see that providing the sharp increase 
of photon Bose condensation critical temperature, the 
argument of logarithm in (57) should tend to unity: 



9a 2 (n - n ph ) 

> 1, T c -> oo, (58) 

that implies the need to fulfill the following relation (see 
d56b): 



fast regime of photons condensation (57). To do it, one 



has to examine the expression ([55]) in the region, where 
no small parameters exist. For this reason we use the nu- 
merical methods as auxiliary. To this end, we introduce 
the following dimensionlcss quantities: 



7 ; 



E a2 — E at — hbjQ 

_ ff *r(3/2) C(3/2) 



.912 



V2ir 2 h?r 



(e a2 — e ai — hu>o)] 



9ai 
9a 2 
3/2 



(60) 



In terms of these dimensionlcss parameters the equation 
(55) can be reduced to the form: 
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|r| 3/2 + (l + .g 12 cxp(l/r))- 1 =7 ? . 



(61) 



Recall again that in accordance with the results, obtained 
in Section [m] the region of negative r is not very promis- 
ing from the viewpoint of BEC implementation at high 
temperatures. Since we have agreed in this section to 
assume (e Q2 — e Ql — hoj ) > 0, the absolute value sign in 



(61 1 can be omitted. 



Figure 2 shows the results of numerical solution of the 
equation (61), reflecting the dependences r (77) for differ- 
ent values of the system parameter 7. These plots, in 
accordance with the definitions (60), reveal the depen- 



dence of the critical temperature of photons Bose con- 
densation on the density of pumped into the system pho- 
tons. To simplify calculations, we assumed g ai /ga 2 = 1- 
The power law dependence (lower solid line) corresponds 
to the numerical solution of equation (61) with 7 = 1, 



i.e. in the case when the first term of equation (55 1 dom 



inates. Logarithmic regime of condensation (upper solid 
line) occurs for values 7 ~ 10 -2 and lower, i.e. in the case 



when the first term in equation (61) (also see (55)) can 



be neglected. Intermediate regimes between "power" and 
"logarithmic", for lucidity, are shown by dashed curves. 

The results of numerical solution witness the transition 
to the fast regime of Bose condensation of photons if 
the parameters of system satisfy the relation 7 < 10~ 2 
(see (60)). Therefore, the following condition is in fact 



imposed on the quantities that characterize the system 
under study 



n p h 



9a\ ~r 9a 2 



(59) 



Note that the expression (59) satisfies the requirement 
n p h ~ n. In other words, the number of photons, pumped 
into the system, is sufficiently large to provide macro- 
scopically large number of free photons in system. Thus, 
in accordance with ( 59 ) the amount of pumped photons 



should be less than the total number of atoms, therefore 
in real experiment the probability of excitation of higher 
atomic levels is extremely small. This circumstance can 
justify the application of the two-level atoms model. 

We study now in details the transition from the tradi- 
tional for Bose condensation "power" regime (51) to the 



nh 3 



[to* (e e 



hio )} 3/2 < 1 



(62) 



In addition, as can be seen from (56), the requirement 
n implies the condition (e a2 — e. 



huj ) 



Thus, the condition for realization the fast condensation 



regime ( 62 ) can be also written as 



T c m* 



3/2 



< 1. 



(63) 



Note that this condition can be achieved in dilute gases 
at rather high temperatures, for example, room temper- 
atures. Thus, for example, for the same photon mass 
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TO* ~ 10 33 g (see [11]), room temperature (T = 300-K"), 
and the density of diluted gas n = 10 16 cm -3 , one obtains 
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-5 



As it is seen, the condition (63) is satisfied 



very well, even up to characteristic temperatures of some 
types of plasma. 



One can easily show that the conditions (34) are also 
satisfied. In fact, there are no doubts in validity of con- 
ditions for atomic components to be nondegenerated, be- 
cause they are valid even in the case of lower tempera- 
tures, determined by power condensation regime (51 ). As 



for the conditions for photonic component to be degener- 
ate, it is implied by the correctness of equality fi* = Itwq 
(see ([37])). 



It is clear from the above analysis that the achievement 
of the "high temperature" regime of Bose condensation 
of photons is essentially defined by both the total den- 
sity of atomic components and the density of the total 
amount of photons, pumped into the system. Thus, one 
can affect the achievement of the transition point both 
with the change of atoms number in system, and with 
pumping photons with a help of laser. Naturally, the 
second option seems much more practical from the point 
of view of experiments. In this regard, the following ques- 
tion rises: what is the rate of achievement the transition 
temperature in the case of some additional pumping of 
photons in the system under study? As a quantitative 
characteristic of the Bose condensation process rate it is 
natural to choose the quantity dT c /dn p h- It should be 
emphasized that we do not actually proceed to the dy- 
namic characteristics of the system. It is assumed that 
after some controlled pumping of additional photons in 
the system (for example, laser impulse pumping) , the sys- 
tem requires some time to come to the thermodynamic 
equilibrium state. However, for the simplicity in subse- 




photons density 



FIG. 2: Dependence of Bose condensation critical tempera- 
ture t on the density of pumped into the system photons r\ 
for different regimes (i.e. for diffe ren t characteristics of sys- 
tem 7), see notations of variables (60 1. For simplicity, all the 
curves are shown for g ai /ga 2 = 1- The lowest solid line cor- 
responds to the true "power" regime ( |51[ ), the highest solid 
line is for the true "logarithmic" regime (57). Intermediate 



regimes for different values of 7 are shown in dashed lines. 



quent considerations, we will call the introduced quantity 
dT c /dn p h the rate of Bose condensation of photons. 

In the case of logarithmic dependence of critical tem- 
perature on the total photon density (see (57)), one can 
readily come to the following expression: 



dT c n (e Q2 - s ai - huj ) r _ 2 



dn 



ph 



n ph (n - n ph ) 



g a2 n - n ph 

9al n ph 



(64) 



In the case, when the "power law" condensation mech- 
anism (given by formula (50)) predominates, one obtains 



dT c 
dn ph 



1 

3m* 



4ir 2 h 3 
T (3/2) C (3/2) 



2/3 



-1/3 -1/3 
n P h xn P h ■ 



(65) 



Comparing expressions ( 64 ) and ( 65 ) it becomes ob- 



viously that the rate of power condensation regime de- 
creases monotonically with increasing the density of 
pumped photons, whereas the rate of logarithmic con- 
densation regime can increase under certain conditions. 
One should simultaneously keep in mind that the require- 
ment n p h ~ n must be held. 

For a more detailed consideration of the introduced 
quantity we will again involve numerical methods. To do 
it, we use the general expression (55) as the equation, 



defining an implicit dependence of the critical temper- 
ature of Bose condensation T c on the total density n p h 
of pumped into the system photons. Taking a derivative 
from the both sides of (55) with respect to n p h, after 
some hackneyed transformations one can obtain: 



d.T 

drj 



7T 1/2 + 



312 exp [1/t] 



r2(l+ 012 exp[l/r]r 



(66) 



Here one should use one more time the definitions ( 60 1 



The obtained expression ( 66 ) should be considered in 
couple with (61) as a system of equations. This sys- 



tem allows to determine both the critical temperature of 
Bose condensation t, and the Bose condensation rate (in 
the sense that was defined above) being characterized by 
quantity dT c /dn p h, as the functions of the pumped pho- 
tons density 77 for given parameters 7 and gi2- 

Figure 3 shows plots for the dependences (77) as a 
result of the numerical solution of equations ( 66 ) and 
(tfjlj) for different values of 7. As in Figure 2, the power 



and logarithmic regimes are shown by the solid lines (bot- 
tom and top ones, respectively), intermediate regimes are 
shown by dashed lines. The plot shows that the rate of 



logarithmic condensation regime ( 64 ) increases with the 



number of photons while the rate of power condensation 
regime (65) decreases. Therefore, the curves in Figure 



3 give an understanding how with an appropriate selec- 
tion of system parameters (see (60)) one can switch the 
system under study to the desired regime with any given 
condensation rate. Note also that at low densities of 



12 




photons density 



FIG. 3: Dependence of condensation rate (r/) on the den- 
sity of pumped into the system pho tons rj for different system 
characteristics 7, see notations (60 1. For simplicity, all curves 
are given for g ai /ga 2 — 1- The lowest solid line corresponds 
to the true "power" regime ( |51[), t he highest solid line is for 
the true "logarithmic" regime~ |57| . Intermediate regimes for 
different values of 7 are sho wn in dashed lines. Minimum 
point for logarithmic regime (641 corresponds to 77 ~ 0.08. 



pumped photons, both dependences (64) - (65) behave 



themselve in a similar way. Both dependences intersect 
each other in the vicinity of the minimum point of ex- 
pression (64) (i.e. approximately for rj ~ 0.08). Thus, 



the numerical results witness in favour of the predomi- 
nant implementation of fast regime of photons Bose con- 
densation in a wide range of pumped into the system 
photons densities. 

In other words, trying to achieve experimentally a state 
with a BEC of photons, that are in thermodynamic equi- 
librium with an ideal gas of effectively two-level atoms 
with the help of set of methods, similar to that described 
in [11, 12], the logarithmic condensation regime occurs 
with the overwhelming probability. Moreover, as it was 
already noted, it applies to a wide range of physical sys- 
tems, that are in equilibrium with photons: from ultra- 
cold gases of alkali metals to some certain types of an 
ideal plasma (for plasma classification types and their 
featured parameters see e.g. [20]). 



VI. CONCLUSION 

The main purpose of the present work, as it has been 
already noted, was to reveal the possibility of BEC exis- 
tence in a gas of photons that are in the state of thermo- 
dynamic equilibrium with the gas of atoms that can re- 
peatedly absorb and reemit photons. Such a mechanism 
of gaining the thermodynamic equilibrium was supposed 
in the present paper to be dominating. The main reason 
for this consideration is the desire to have some resem- 
blance between the system under study and the system 
of [11, 12], where photon BEC was first obtained experi- 
mentally. We managed to demonstrate the possibility of 



achieving such a state of the photon component within 
a framework of fairly simple model, indicating, however, 
the conditions under which the situation could be close 
to real one. In some cases, our model requires improve- 
ments or corrections. It is connected, for example, with 
taking into account that Bose condensation of photons 
is affected by interactions between the atoms of atomic 
components and photon scattering by atoms. This is a 
separate, sophisticated issue that in our view is beyond 
the scope of the present paper. Its solution has to be 
in framework of the microscopic approach based on the 
low-temperature quantum electrodynamics. In our view, 
a very promising is the approach proposed in [21]. Au- 
thors are currently working on these issues. 

The atoms of atomic components in the present paper 
are considered to be two-level. It is a fairly widespread 
technique in theoretical physics, optics, and photonics. 
But it should be noted that from the standpoint of this 
paper, it is not a fundamental limitation. It only sim- 
plifies calculations, allowing in some cases to obtain the 
results in analytical form. The equations (15) -(17) can 



be formulated for an arbitrarily large number of com- 
ponents. But in this case there rises a challenge of a 
controlled numerical solution of such equations. 

It is also worth to discuss the problem of photon dis- 
persion in matter. The importance of solving this prob- 
lem is related, first of all, to the fact that the explicit 
form of ui (k) determines the effective mass of a photon 
in the matter (see (7) - (13)). In this paper we have used 
for numerical estimations the value of the effective mass 
of the photon from [11]. However, the effective mass of 
photons in different substances can significantly vary. It 
can be related to the interaction features of photons with 
the structural units of a particular substance. For this 
reason, ideally the dependence of the photon frequency 
on the wave vector uj (k) has to be calculated indepen- 
dently for each specific system in a self-consistent way. 
Moreover, such self-consistent calculations of the photon 
dispersion in matter must proceed from the first princi- 
ples underlying the description of quantum many-particle 
systems. In other words, in the case of issues related to 
the description of BEC of photons, calculations of photon 
dispersion in matter have to be based on nonrelativistic 
quantum electrodynamics and quantum statistical me- 
chanics. In its most general form, this problem seems to 
be unreasonable. However, in some cases, it is quite re- 
alistic. It was successfully accomplished, for example, in 
[17,18] for the description of slowing and absorption pro- 
cesses of electromagnetic waves in gases of alkali metals 
with BEC or propagation of relativistic charged parti- 
cles in the same substances [19]. These studies used a 
microscopic approach based on the low energy quantum 
electrodynamics, constructed in [21]. There is a hope to 
generalize the approach of [21], [17-19] in order to cal- 
culate the effective mass of the photon in the systems 
similar to those considered in the present paper. 

In conclusion, the authors would also emphasize that 
mentioning an ideal plasma in our work is related to the 
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demonstration of the vastness of the temperature range these results to photons in plasma. For a description of 
in which used approach and obtained results are valid, such processes in an ideal plasma, our model requires a 
but it is not a statement about the direct applicability of substantial modification (in this regard see [9]). 



[1] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wie- 

man, E.A. Cornell, Science 269 (1995), 198. 
[2] K.B. Davis, M.-O. Mewes, M.R. Andrews, N.J. van 

Druten, D.S. Durfee, D.M. Kurn, W. Ketterle, Phys. 

Rev. Lett. 75, (1995), 3969. 
[3] C.C. Bradley, C.A. Sackett, J.J. Tollett, R.G. Hulet, 

Phys. Rev. Lett., 75, (1995), 1687. 
[4] E.A. Donley, N.R. Claussen, S.T. Thompson, and C.E. 

Wieman, Nature 417 (2002), 529. 461. 
[5] S. Durr, T. Volz, A. Marte, and G. Rempe, Phys. Rev. 

Lett. 92, (2004), 020406. 
[6] A. Einstein, Sitzungsber. Kgl. Preuss. Akad. Wiss., 1 

(1925), 3. 

[7] S.N. Bose, Z. Phys., 26 (1924) 178. 

[8] V.A. Kuz'min, M.E. Shaposhnikov, JETP Letters, 27 
(1978) 628. 

[9] Y.B. Zel'dovich, E.V. Levich, Sov. Phys. JETP, 28 (1969) 
1287. 

[10] R.Y. Chiao. and J. Boyce, Phys. Rev., A 60, (1999) 4114. 
[11] J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Nature, 

468 (2010) 545. 
[12] J. Klaers, J. Schmitt, T. Damm, D. Dung, F. Vewinger, 



and M. Weitz, arXiv:1303.5772vl [cond- mat. quant-gas]. 
[13] L.D. Landau, E.M. Lifshitz, Statistical Physics, Part 1. 

Vol. 5 (2nd ed.). ((Pergamon, Oxford, 1970)) 
[14] Yu.V. Slyusarenko,and A.G. Sotnikov, Low Temp. Phys., 

33, 30 (2007). 

[15] A.I. Akhiezer, S.V. Peletminskii, Methods of Statistical 

Physics (Pergamon, Oxford, 1981). 
[16] Y. V. Slyusarenko and A. G. Sotnikov, Condens. Matter 

Phys. 9, 459 (2006). 
[17] Y. Slyusarenko and A. Sotnikov, Phys. Rev. A 78, 053622 

(2008) . 

[18] Y. Slyusarenko and A. Sotnikov, Phys. Rev. A 80, 053604 

(2009) . 

[19] Y. Slyusarenko and A. Sotnikov, Phys. Rev. A 83, 023601 
(2011). 

[20] Peter A. Sturrock, Plasma Physics: An Introduction to 
the Theory of Astrophysical, Geophysical and Labora- 
tory Plasmas (Cambridge University Press, 1994). 

[21] S.V. Peletminskii, Yu.V Slyusarenko, J. Math. Phys., 46, 
022301 (2005). 



